Two-proton radioactivity and three-body decay. III. Integral formulae for decay 
widths in a simplified semianalytical approach. 



L. V. Grigorenko 1 - 2 - 3 and M. V. Zhukov 4 

1 Flerov Laboratory of Nuclear Reactions, JINR, RU-14-1980 Dubna, Russia 
2 Gesellschaft fur Schwerionenforschung mbH, Planckstrasse 1, D-64291, Darmstadt, Germany 
1 1 RRC "The Kurchatov Institute", Kurchatov sq. 1, 123182 Moscow, Russia 
^Fundamental Physics, Chalmers University of Technology, S-^1296 Goteborg, Sweden 

Three-body decays of resonant states are studied using integral formulae for decay widths. Theo- 
retical approach with a simplified Hamiltonian allows semianalytical treatment of the problem. The 
model is applied to decays of the first excited 3/2~ state of 17 Ne and the 3/2~ ground state of 45 Fe. 
The convergence of three-body hyperspherical model calculations to the exact result for widths and 
energy distributions are studied. The theoretical results for 17 Ne and 45 Fe decays are updated and 
uncertainties of the derived values are discussed in detail. Correlations for the decay of 17 Ne 3/2~ 
state are also studied. 
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emission, 21.10.Tg — Lifetimes 



I. INTRODUCTION 



(a) 



(b) 



The idea of the "true" two-proton radioactivity was 
proposed about 50 years ago in a classical paper of 
Goldansky The word "true" denotes here that we 
are dealing not with a relatively simple emission of two 
protons, which becomes possible in every nucleus above 
two-proton decay threshold, but with a specific situa- 
tion where one-proton emission is energetically (due to 
the proton separation energy in the daughter system) or 
dynamically (due to various reasons) prohibited. Only 
simultaneous emission of two protons is possible in that 
case (see Fig. [TJ more details on the modes of the three- 
body decays can be found in Ref. (H). The dynamics of 
such decays can not be reduced to a sequence of two-body 
decays and from theoretical point of view we have to deal 
with a three-body Coulomb problem in the continuum, 
which is known to be very complicated. 

Progress in this field was quite slow. Only recently a 
consistent quantum mechanical theory of the process was 
developed @, H, 0] , which allows to study the two-proton 
(three-body) decay phenomenon in a three-body cluster 
model. It has been applied to a range of a light nuclear 
systems ( 12 0, 16 Ne g, 6 Be, 8 Li*, 9 Be* @, 17 Ne*, 19 Mg 
[7|). Systematic exploratory studies of heavier prospec- 
tive 2p emitters 30 Ar, 34 Ca, 45 Fe, 48 Ni, 54 Zn, 58 Ge, 62 Se, 
and 66 Kr 0, Q) have been performed providing predic- 
tions of lifetime ranges and possible correlations among 
fragments. 

Experimental studies of the two-proton radioactivity is 
presently an actively developing field. Since the first ex- 
perimental identification of 2p radioactivity in 45 Fe [9l[l0j 
it was also found in 54 Zn [1JJ. Some fingerprints of the 
48 Ni 2p decay were observed and the 45 Fe lifetime and de- 
cay energy were measured with improved accuracy 
There was an intriguing discovery of the extreme en- 
hancement of the 2p decay mode for the high-spin 21 + 
isomer of 94 Ag, interpreted so far only in terms of the 
hyperdeformation of this state [l3j . New experiments, 




Three-body g 
decay 
window 



A (A-2) + 2N A (A-2) + 2N 

(A-l) + N thresholds (A-l) + N 



FIG. 1: Energy conditions for different modes of the two- 
nucleon emission (three-body decay): true three-body decay 
(a), sequential decay (b). 



aimed at more detailed 2p decay studies (e.g. observa- 
tion of correlations), are under way at GSI ( 19 Mg), MSU 
( 45 Fe), GANIL ( 45 Fe), and Jyvaskyla ( 94 Ag). 

Several other theoretical approaches were applied to 
the problem in the recent years. We should mention 
the "diproton" model [3, [HI, "R-matrix" approach 
[l|| [13, El, HH, continuum shell model [2(|, and adiabatic 
hyperspherical approach of [21]. Some issues of a com- 
patibility between different approaches will be addressed 
in this work. 

Another, possibly very important, field of applica- 
tion of the two-proton decay studies was shown in Refs. 
[22l [23| . It was demonstrated in [12] that the importance 
of direct resonant two-proton radiative capture processes 
was underestimated in earlier treatment of the rp- process 
waiting points [H| . The scale of modification of the astro- 
physical 2p capture rates can be as large as several orders 
of magnitude in certain temperature ranges. In paper (23j 
it has been found that nonresonant El contributions to 
three-body (two-proton) capture rates can also be much 
larger than was expected before. The updated 2p as- 
trophysical capture rate for the 15 0(2p,7) 17 Ne reaction 
appears to be competing with the standard 15 0(a,7) 19 Ne 
breakout reaction for the hot CNO cycle. The improve- 
ments of the 2p capture rates obtained in [22|, HH are 
connected to consistent quantum mechanical treatment 
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of the three-body Coulomb continuum in contrast to the 
essentially quasiclassical approach typically used in as- 
trophysical calculations of three-body capture reactions 
(e.g. MM)- 

The growing quality of the experimental studies of the 
2p decays and the high precision required for certain as- 
trophysical calculations inspired us to revisit the issues 
connected with different uncertainties and technical diffi- 
culties of our studies. In this work we make the following, 
(i) Extend the two-body formalism of the integral for- 
mulae for width to the three-body case. We perform the 
relevant derivations for the two-body case to make the 
relevant approximations and assumptions explicit, (ii) 
Formulate a simplified three-body model which has many 
dynamical features similar to the realistic case, but allows 
the exact semianalytical treatment and thus makes pos- 
sible a precise calibration of three-body calculations, ft 
is also possible to study in great detail several impor- 
tant dependencies of three-body widths in the frame of 
this model, (iii) Perform practical studies of some sys- 
tems of interest and demonstrate a connection between 
the simplified semianalytical formalism and the realistic 
three-body calculations. 

The unit system fi — c — 1 is used in the article. 



where v = \/2E/M. The pole solution for Hamiltonian 



(H - E r )^+\r) = (T + V - E r )^(r) = ° 
provides the WF with outgoing asymptotic 
^(r)=r- 1 ^ + \kr)Y lm (r). 
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(1) 



For single channel two-body problem the pole solution is 
formed only for one selected value of angular momentum 
I. In the asymptotic region 

(k r r) r = H H[ +) (k r r) = Gi (k r r) + iF t (k r r) . (2) 
The above asymptotic is growing exponentially 

ipi + \k r r) r ~ fl exp[+ifc r r] w exp[+ifc r r] exp[+rr/(2tv)] 

as a function of the radius at pole energy. This unphys- 
ical growth is connected to the use of time-independent 
formalism and could be reliably neglected for typical ra- 
dioactivity time scale as it has a noticeable effect at very 
large distances. 

Applying Green's procedure to complex energy WF 



*(+)t 



(H - £? r )*(+)l - \(H - E r )¥ + )] = 



II. INTEGRAL FORMULA FOR WIDTH 

Integral formalisms of width calculations for narrow 
two-body states are known for a long time, e.g. [26l. [27|] . 
The prime objective of those studies was a-decay widths. 
An interesting overview of this field can be found in the 
book [281 ] . This approach, to our opinion, did not pro- 
duce novel results as the inherent uncertainties of the 
method are essentially the same as those of the R-matrix 
phenomenology, which is technically much simpler (see 
e.g. a discussion in [29]). An important nontrivial appli- 
cation of the integral formalism was calculation of widths 
for proton emission off deformed states [3(| HH • There 
were attempts to extend the integral formalism to the 
three-body decays, using a formal generalization for the 
hyperspherical space U [32j. These were shown to be 
difficult with respect to technical realisation and to be 
inferior to other methods developed in [1, Q . 

Here we develop an integral formalism for the three- 
body (two-proton) decay width in a different way. How- 
ever, first we review the standard formalism to define 
(clearer) the approximations used. 



we get for the partial components at pole energy E r 



zr^ (+) *^ (+) 
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2M 



air 2 



^2, /,( + )* 

d y>i ,/.(+) 



After radial integration from to R (here and below R 
denotes the radius sufficiently large that the nuclear in- 
teraction disappears) we obtain 



(+) 



2Mz f R 



(+) 



dr 



2L 



which corresponds to a definition of the width as a decay 
probability (reciprocal of the lifetime): 

N = N cxp[-</r] = N Q exp[~rt] . 

The width T is then equal to the outgoing flux ji through 
the sphere of sufficiently large radius R, divided by num- 
ber of particles Ni inside the sphere. 

Using Eq. (f5|) the flux in the asymptotic region could 
be rewritten for k r — » k r in terms of a Wronskian 



A. Width definition, complex energy WF 

For decay studies we consider the wave function (WF) 
with complex pole energy 

E r = ~k 2 r /{2M) = E r - iT/2 , ~k r 7ak r - iT/(2v r ) , 



1 



Ji 



2Mi 



(+)« 



dr 



(+) 



Tr^ 



(+)» 



(+) 



r=R 

(4) 



= (k r /M) W(Fi(k r R), Gi(k r R)) = v r , 
where the Wronskian for real energy functions F, \ , Gi is 
W(F l .,G l ) = G l F' l -G' l F l = l . 
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The effect of the complex energy is easy to estimate (ac- 
tually without loss of a generality) in a small energy ap- 
proximation 



kr^O „ i + 1 „ fer— (kr) 1 



F t (kr) ~ C l (krY + \ G t (kr) 



(21 + 1)C, 



(5) 



where C; is a Coulomb coefficient (defined e.g. in Ref. 
The flux is then 
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k r H\ \k*r)H[ +) \k r r)-k;H\- ] \k* r r)H\ + \k r r) 
2iM 



So, the equality ([4]) is always valid for I — and for / ^ 

we get 



i(J + i) 



1/2 



E r 



B. Two-body case, real energy WF 

Now we need a WF as real energy E — k 2 /2M solution 
of Schrodinger equation 

(H - J5)* k (r) = (T + y nuc + V coul - #)# k (r) = , 
* k (r)=47rV i'(Ar)-Vi(*r)X; 5ft,(*)*W*) , 

* — *i * — 

in S-matrix representation, which means that for r > R 

Mkr) = % - [(Gi(kr) - iF t (kr)) - 5i(G,(*t) + »F,(fcr))] . 

At resonance energy £? r 

S,(£V) = e 2l5i ^ = e 2 "/ 2 = -1 

and in asymptotic region, defined by the maximal size of 
nuclear interaction R, 



ipi(k r r) r = fl i Gi(k r r) 



At resonance energy we can define a "quasibound" WF 
ipi as matching the irregular solution Gi and normalized 
to unity for the integration in the internal region limited 
by radius R: 

^ (M - — — = -t — — . (6) 



f/.i: 



iV, 



1/2 



Now we introduce an auxiliary Hamiltonian H with 
different short range nuclear interaction V nuc , 



and also construct other WF in S-matrix representation 

* — * / * — ' m 

W (fcr) = ~[(Gi(kr)-iF l (kr))-Si(G l (kr)+iF l (kr))] 

for r > R. Or in equivalent form: 

¥ >,(*T)=exp(«,) [F(fcr)cos(5 / ) + G / (A;r)sm(^)] . (7) 

The Hamiltonian H should provide the WF $k(r) which 
at energy E r is sufficiently far from being a resonance 
WF and for this WF S t (E r ) ~ 0. 

For real energy WFs \&k(r) and $k(r) we can write: 



* k (r)t [(H - E)* k (r)] - [{H - E)*k(r)] f * k (r) = . 



,(8) 



For WFs taken at resonance energy E r this expression 
provides 

i-R pR. 

2M tf{V -V)ipidr = 2MiN t 1/2 <p\ (V - V)ipidr 
Jo Jo 

= exp(-*J,) oos(«,) K W(F l (k r R),Gi(k r R)) ,(9) 
-i exp(-iSi) cos((5;) fc r 



A7 



1/2 



2M f R ^(y-^Mdr 



From Eqs. Q, (f4]), ([6|) and the approximation 
it follows that 



(+) 



r 



(+) 



'dr Jo Ml dr 



1/2 



iV cos 2 (<5/) 



^(V-y)^dr 



(10) 



So, the idea of the integral method is to define the in- 
ternal normalizations for the WF with resonant boundary 
conditions (this is equivalent to determination of the out- 
going flux for normalized "quasibound" WF) by the help 
of the eigenfunction of the auxiliary Hamiltonian, which 
has the same long-range behaviour and differs only in the 
compact region. 



III. ALTERNATIVE DERIVATION 

Let us reformulate the derivation of Eq. (JTTJJ) in a more 
general way, so that the detailed knowledge of the WF 
structure for ipi and ^. is n °t required. It would 
allow a straightforward extension of the formalism to 
the three-body case. We start from Schrodinger equa- 
tion in continuum with solution at the pole energy 
E r = E r +ir/2: 



(H - E)$ k (r) = (T + V nuc + V coul - £7)* k (r) = , (h - E r ) = (t + V - E r ) = . (11) 
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Then we rewrite it identically via the auxiliary Hamilto- 
nian H = T + V 

(h + v-v-e^) * (+) = {v -V) * (+) 

(H~E r )¥ +) = (V -V + iT/2)^> {+) . (12) 

Thus we can use the real-energy Green's function Ge t of 
auxiliary Hamiltonian H to "regenerate" the WF with 
outgoing asymptotic 



= GeJ (V-V + i r/2) $ 



(+) 



(13) 



At this point in Eq. (JT3J) * (+) _= * (+) and the bar in the 
notation for "corrected" WF $(+) is introduced for later 
use to distinguish it from the "initial" WF [the one 
before application of Eq. IJT3J)]. Further assumptions we 
should consider separately in two-body and three-body 
cases. 



A. Two-body case 

To define the width T by Eq. © we need to know 
the complex-energy solution at pole energy. For 

narrow states r <C E r this solution can be obtained in a 
simplified way using the following approximations. 

(i) For narrow states we can always choose the auxil- 
iary Hamiltonian in such a way that r <C V — V , and we 
can assume r — > in the Eq. (fl~3]) . 

(ii) Instead of complex-energy solution $(+) in the 
right-hand side of (fToT) we can use the normalized real- 
energy quasibound solution vff defined for one real reso- 
nant value of energy E r = k 2 /2M 



N, 
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j dQ, J dr r 2 


*im(r) 


= 1 









So, the Eq. (fl"3j) is used in the form 



(+) 

lm 



The solution \I>(+) is matched to function 
h\ + \kr) = G t (kr) + iFi(kr) 



(14) 
(15) 



while the solution ^> is matched to function Gi . For deep 
subbarrier energies it is reasonable to expect that in the 
internal region r < R 



Gi^>Fi 



Re[* 



(+)■ 



* > Im[# 



/(+)1 



In the single channel case it can be shown by direct cal- 
culation that an approximate equality 



MR 2 T 



Re[# 



(+)• 



> 



Im[# 



(+)1 



holds in the internal region and thus for narrow states 
r <C E r the approximation — > ([T4")) should be very 
reliable. 
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FIG. 2: Single particle coordinate systems: (a) "V" system 
typical for a shell model. In the Jacobi "T" system (b), 
"diproton" and core are explicitly in configurations with def- 
inite angular momenta l x and l y . For a heavy core the Jacobi 
"Y" system (c) is close to the single particle system (a). 



To derive Eq. (fTU)) the WF with outgoing asymptotic 
is generated using the Green's function of the auxiliary 
Hamiltonian H and the "transition potential" (V — V). 
The standard two-body Green's function is 



G 



(+) 

k 2 /(2r. 



,(r,r') 



2My-J ipi(kr)h\ +) (kr'),r < r ' 1 
krr' 2f \ h [ , + \kr)ipi{kr'),r> r' J 

x Y, Y lm (r)YL(r) , (16) 

where the radial WFs h\ and tpi of the auxiliary Hamil- 
tonian are defined in (|15p and ([7]). 



*, (+) (r) 

lm V / 



dr 



'G^ /{2m) {v,v'){V-V) ^w(r') 



For the asymptotic region r > R 

= ^^ (+) (*rr)H m (f) 
ii 



x / dr' ^(fc r rO (F - V") tpi(kr,r') 
Jo 

The outgoing flux is then calculated [see Eq. (j4|)] 



31 = ^—Y I dn — lm *, (+) *(r)V*, ( " 
21 + 1^ M I lm w ir ' 



r=_R 



As far as function \I/ is normalized by construction then 

2 



r = ji = - 



R 







dr<pi(k r r) (V~V)Mk r ,r) 



(17) 



Note, that this equation differs from Eq. (|10[) only by a 
factor l/(cos 2 [<5/]) which should be very close to unity for 
sufficiently high barriers. 



B. Simplified model for three-body case 

In papers 0, Q the widths for three-body decays were 
defined by the following procedure. We solve numerically 
the problem 



{H-E 3r )V = 
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with some box boundary conditions (e.g. zero or qua- 
sibound in diagonal channels at large distances) getting 
the WF ^> normalized in the finite domain and the value 
of the real resonant energy E 3r . Thereupon we search for 
the outgoing solution of the equation 

with approximate boundary conditions of three-body 
Coulomb problem (see Ref. [2| for details) and arbitrary 
r. The width is then defined as the flux through the 
hypersphere of the large radius divided by normalization 
within this radius: 



f<m 5 ^(+)* p 5/2jL (0 5/2^(+) 



r = = 

N 



P=Pn 



m J <m 5 J pmax P 5 dp |*(+)| 



The 3-body WF with outgoing asymptotic is 



(18) 



(19) 



where the definitions of the hyperspherical variables p, 

J.JM 
>K~i 



O5 and hyperspherical harmonics Ji^ can be found in 



Rcf. g. 

Here we formulate the simplified three-body model in 
the way which, on one hand, keeps the important dy- 
namical features of the three-body decays (typical sizes 
of the nuclear potentials, typical energies in the subsys- 
tems, correct ratios of masses, etc.), and, on the other 
hand, allows a semianalytical treatment of the problem. 
Two types of approximations are made here. 

The three-body Coulomb interaction is 



ycoul 



X 



Z\Z 3 OL 



Z 2 Z 3 a 



Y^ 


A 2 X 


- 1 " 


V MX. 




A t +A 2 




1 A t +A 2 



(20) 



where a is the fine structure constant. By convention, see 
e.g. Fig. [2J in the "T" Jacobi system the core is particle 
number 3 and in "Y" system it is particle number 2. We 
assume that the above potential can be approximated by 
Coulomb terms which depend on Jacobi variables X and 
Y only: 



Z x a 
~X~ 



T/-COU] 

y 



(Y) 



Z y a 



(in reality for the small X and Y values the Coulomb 
formfactors of the homogeneously charged sphere with 
radius r sp h are always used). The effective charges Z x 
and Z y could be considered in two ways. 

1. We can neglect one of the Coulomb interactions. 
This approximation is consistent with physical situ- 
ation of heavy core and treatment of two final state 
interactions. Such a situation presumes that Jacobi 
"Y" system is preferable and there is a symmetry 
in the treatment of the X and Y coordinates, which 
are close to shell-model single particle coordinates. 



Z\Z C . 



Zr>Zr 



(21) 



Further we refer this approximation as "no p- 
p Coulomb" case, as typically the proton-proton 
Coulomb interaction is neglected compared to 
Coulomb interaction of a proton with heavy core. 

2. We can also consider two particles on the X coor- 
dinate as one single particle. The Coulomb inter- 
action in p-p channel is thus somehow taken into 
account effectively via a modification of the Z y 
charge: 



Z x — Z\ Z, 



1 ^core 



Zy — Zl{Z c 



Zx) 



(22) 



Below we call this situation as "effective p-p 
Coulomb" case. 

For nuclear interactions we can assume that 

1. There is only one nuclear pairwise interaction and 

H = T + V 3 (p) + V™ u \X) 
+V™ C {X) + V y cou \Y) , 
AV(X,Y) = V™ C (Y)-V 3 (p). (23) 



This approximation is good for methodological pur- 
poses as it allows to focus on one degree of freedom 
and isolate it from the others. From physical point 
of view it could be reasonable if only one FSI is 
strong [42j], or we have reasons to think that de- 
cay mechanism associated with this particular FSI 
is dominating. Potential V? UC (Y) in the auxiliary 
Hamiltonian (|27| is "unphysical" in that case and 
can be put zero [43|]. We further refer this model 
as "one final state interaction" (OFSI). 

2. We can consider two final state interactions (TFSI). 
Simple form of the Green's function in that case can 
be preserved only if the core mass is considered as 
infinite (the X and Y coordinates in the Jacobi 
"Y" system coincide with single-particle core-p co- 
ordinates). In that case both pairwise interactions 
T/ nuc (Y) and V^ nuc (Y) are treated as "physical", 
that means that they are both present in the ini- 
tial and in the auxiliary Hamiltonians. Thus only 
three-body potential "survive" the V — V subtrac- 
tion: 



H = T + V 3 (p) + V a 



Hx) + v x ™ c (x) 

\Y) + V™ C (Y) , 
AV(X,Y) = -V 3 (p) . 



+V y cou 



(24) 



The three-body potential is used in this work in 
Woods-Saxon form 



V 3 (p) = V 3 °(l + c W [(p-p )/a p }y 



(25) 



with po — 5 fm for 17 Ne, po = 6 fm for 45 Fe [Hj], and 
a small value of diffuseness parameter a p — 0.4 fm. Use 
of such three-body potential is an important difference 
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from our previous calculations, where it was utilized in 
the form 



v 3 (p) = v°(i + (p/p ) 3 y 



(26) 



which provides the long-range behaviour ~ P ■ Such 
an asymptotic in p variable is produced by short-range 
pairwise nuclear interactions and thus the interpretation 
of three-body potential (|2"6"|) is phenomenological taking 
into account those components of pairwise interactions 
which were omitted for some reasons in calculations. In 
this work the aim of the potential V3 is different. On 
one hand we would like to keep the three-body energy 
fixed while the properties (and number) of pairwise in- 
teractions are varied. On the other hand we do not want 
to change the properties of the Coulomb barriers beyond 
the typical nuclear distance (this is achieved by the small 
diffuscness of the potential) . Thus this potential is phe- 
nomenological taking into account interactions that act 
only when both valence nucleons are close to the core 
(both move in the mean field of the nucleus). 

The auxiliary Hamiltonian is taken in the form that 
allows a separate treatment of X and Y variables 



H 



T + V/ coul (X ) + V™ c (X ) + V™ ul (Y ) + V y nuc (Y) (27) 



In this formulation of the model the Coulomb potentials 
are fixed as shown above. The nuclear potential V^ nuc (X) 
|V y nuc (F) if present] defines the position of the state in the 
X [Y] subsystem. The three-body potential V 3 (p) defines 
the position of the three-body state^which is found using 
the three-body HH approach of (2, 4]. After that a new 
WF with outgoing asymptotic is generated by means of 
the three-body Green's function which can be written for 
(1271) in a factorized form (without paying attention to the 
angular coupling) 

G^l (XY, X'Y') = — J dE x Gj+> (X, X') G« (Y, Y'), 

where E 3r = E x + E y [E x , E x are energies of subsys- 
tems). The two-body Green's functions in the expres- 
sions above are defined as in (fTo| via eigenfunctions of 
the subhamiltonians 

S X -E X =T X + V x coul (X) + V™ C (X) - E x 

Hy - Ey =Ty + (Y) + V™ (Y) - Ey ' 

In the OFSI case the nuclear potential in the "Y" sub- 
system should be put V; nuc (Y) = 0. The "corrected" 
continuum WF \ff( + ) is 
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FIG. 3: Convergence of the 17 Ne width in a simplified model 
in the "Y" Jacobi system. One final state interaction model 
with experimental position E-z r — 0.535 KeV of the s-wave 
two-body resonance. Diamonds show the results of dynamic 
HH calculations. Solid curves correspond to calculations with 
effective FR potentials. 



The asymptotic form of the "corrected" continuum 
WF is 



V ./A/(X Y ) - — 



de 



A(e) 



2-Ki XY J v x {e)v y (e) 

Xe ik x (e)X+ik v (e)Y ^ g> y £ ,g £] ^ 

E x = eE 3r ; E y = (1 - e)E 3r ; Vl (e) = yj2E l jM l 
A(e) = [ dX' [ dY' ^(k x {e)X') Vly (k y {e)Y') 



xAV(X',Y') <pu.i v s(X',Y').{29) 

The "corrected" outgoing flux j c can be calculated on the 
sphere of the large radius for any of two Jacobi variables. 
E.g. for X coordinate we have [45| 



j c (E 3r ) = Im 



X 2 I dfij dY(*W*^f(+) 



E 3r 



2ir v x v y Jo M x 2ir 



x 2 7rS(k y (e')-k y (e)) 



(30) 



Values v[ above denote Vi(s'). The flux is obtained as 



2ni 

■,(+) 



*W(X,Y) - / dX'dY' I dE x G { +\x,X r ) 



x G^(Y,Y') AV(X',Y') f (+) (X'Y') 

The "initial" solution ^ (+) of Eq. (19]) rewritten in the 
coordinates X and Y is 



<!>(+) (1 r ~Y\ - ^LUySjX, Y) 

w jmI^' Y i 



[[ly®Q L ® S ]jM ( 28 ) 



j c {E 3r ) = -E 3r de \A(e)\ 2 . (31) 

7T Jo v x {e)v y (e) 

In principle as we have seen above that the widths ob- 
tained with both fluxes Eqs. (fT5)l and (f3"Tj) should be 
equal 
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FIG. 4: Convergence of widths in OFSI model for different 
positions Ei r of the two-body resonance in the core-p channel 
(Jacobi "Y" system). For Jf max > 24 the value of if max de- 
note the size of the basis for Feshbach reduction to 7-Gnax = 24. 

This is the idea of calibration procedure for the simplified 
three-body model. The convergence of the HH method 
(for WF ^j^-) is expected to be fast in the internal re- 
gion and much slower in the distant subbarrier region. 
This should be true for the width T calculated in the HH 
method. However, the procedure for calculation of the 
"corrected" width r c is exact under the barrier and it is 
sensitive only to HH convergence in the internal region, 
which is achieved easily. Below we demonstrate this in 
particular calculations. 

IV. DECAYS OF THE 17 NE 3/2 AND 45 FE 3/2~ 
STATES IN A SIMPLIFIED MODEL 

In this Section when we refer widths of 17 Ne and 45 Fe 
we always mean the 17 Ne 3/2" state (E 3r = 0.344 MeV) 
and the 45 Fe 3/2" ground state (E 3r = 1.154 MeV) cal- 
culated in a very simple models. We expect that impor- 
tant regularities found for these models should be true 
also in realistic calculations. However, particular values 
obtained in realistic models may differ significantly, and 
this issue is considered specially in the Section fVl 

To keep only the most significant features of the sys- 
tems we assume pure sd structure (l x — 0, l y = 2) for 
17 Ne and pure p 2 structure (l x = l,l y — 1) for 45 Fe in 
"Y" Jacobi system (see Fig. [2]). Spin dependencies of 
the interactions are neglected. The Gaussian formfactor 

Vr c (r)=V i0 e X p[-(r/ro) 2 } , 

where i = {x,y}, is taken for 17 Ne (see Table HJ, and a 
standard Woods-Saxon formfactor is used for 45 Fe (see 
Table 10, 

Vr c (r) = V i0 [1 + exp[(r - r^/a]}- 1 . (33) 

The simplistic structure models can be expected to 
overestimate the widths. There should be a considerable 



FIG. 5: Width of the 17 Ne 3/2" state as a function of two- 
body resonance position E% r - Dashed, dotted and solid lines 
show cases of pure HH calculations with Jfmax = 24, the 
same but with Feshbach reduction from if max = 100, and the 
corrected width F c . Inset shows the same, but as a function 
of the potential depth parameter V x o- Gray area shows the 
transition region from three-body to two-body decay regime. 
The gray curve shows simple analytical dependence of Eq. 
(El- 
weight of d 2 component (l x = 2, l y = 2) in 17 Ne and f 2 
component (l x = 3, l y = 3) in 45 Fe. Also the spin-angular 
coupling should lead to splitting of the single-particle 
strength and corresponding reduction of the width es- 
timates (e.g. we assume one s-wave state at 0.535 keV in 
the "X" subsystem of 17 Ne while in reality there are two 
s-wave states in 16 F: 0" at 0.535 MeV and 1" at 0.728 
MeV). Thus the results of the simplified model should 
most likely be regarded as upper limits for widths. 

A. One final state interaction — core-p channel 

First we take into account only the 0.535 MeV s-wave 
two-body resonance in the 16 F subsystem (this is the 
experimental energy of the first state in 16 F). Conver- 
gence of the 17 Ne width in a simplified model for Jacobi 
"Y" system is shown in Fig. [3] The convergence of the 
corrected width r c as a function of if ma x is very fast: 
K max > 8 for the width is stable within ~ 1%. For maxi- 
mal achieved in the fully dynamic calculation K max = 24 
the three-body width T is calculated within 30% preci- 
sion. Further increase of the effective basis size is possi- 
ble within the adiabatic procedure based on the so called 
Feschbach reduction (FR). 

Feschbach reduction is a procedure, which eliminates 
from the total WF f = f p + $ g an arbitrary subspace q 
using the Green's function of this subspace: 

Hp = T p + V p + VpqGqVpq 

In a certain adiabatic approximation we can assume that 
the radial part of kinetic energy is small and constant un- 
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FIG. 6: Convergence of energy distribution for Ne in the 
"Y" Jacobi system. 



der the centrifugal barrier in the channels with so high 
centrifugal barrier that it is much higher than any other 
interaction. In this approximation the reduction proce- 
dure becomes trivial as it is reduced to construction of 
effective three-body interactions V^L K <„, by matrix op- 



TABLE I: Parameters for 17 Ne calculations. Potential pa- 
rameters for 15 0+p channel in s-wave (V^o in MeV, tq — 3.53 
fm) and 16 F+p channel in a!- wave (V y o in MeV). Radius of the 
charged sphere is r sp h — 3.904 fm. Widths r 4 of the state in 
the subsystem and experimental width values r exp for really 
existing at these energies states are given in keV. The cor- 
rected three-body width r c is given in the units 10 -14 MeV. 
TFSI calculations with d-wave state at 1.2 MeV are made 
with s-wave state at 0.728 MeV. 
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FIG. 7: Energy distributions for 17 Ne in the "Y" Jacobi 
system for different two-body resonance positions E% T . The 
three-body decay energy is E^ r = 0.344 MeV. The distri- 
butions are normalized to have unity value on maximum of 
three-body components. The values near the peaks show the 
fraction of the total intensity concentrated within the peak. 
Note the change of the scale at vertical axis. 
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(if + 3/2) (if + 5/2) 



Et — E 
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Summation over indexes with bar is made for eliminated 
channels. No strong sensitivity to the exact value of the 
"Feshbach energy" Et is found and we take it as Ef = E 
in our calculations. More detailed account of the pro- 
cedure applied within HH method can be found in Ref. 

It can be seen in Fig. [3] (solid line) that Feschbach re- 
duction procedure drastically improves the convergence. 
However, the calculation converges to a width value, 
which is somewhat smaller than the corrected width value 
(that should be exact). The reason for this effect can 
be understood if we make a reduction to a smaller "dy- 
namic" basis size (K max — 12, gray line). The calculation 
in this case also converges, but even to a smaller width 
value. We can conclude that FR procedure allows any- 
how to approach the real width value, but provides a 
good result only for sufficiently large size of the dynamic 
sector of the basis. 

The next issue to be discussed is a convergence of the 
width in calculations with different positions i?2r of two- 
body resonance in the core+p subsystem. It is demon- 
strated for several energies E^ r in Fig. 2J When the 
resonance in the subsystem is absent (or located rela- 
tively high) the convergence of the width value to the 
exact result is very fast both in the pure three-body 
and in the "corrected" calculation (in that case, how- 
ever, much faster). Here even FR is not required as the 
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FIG. 8: Convergence of the 17 Ne width in a simplified model. 
Jacobi "Y" system. OFSI model with s-wave two-body reso- 
nance at i?2r = 0.535 MeV; Gaussian potential and potential 
with repulsive core. TFSI model with d-wave two-body reso- 
nance at E 2r = 0.96 MeV. 



convergent result is achieved in the HH calculations by 
K max = 10 — 24. The closer two-body resonance ap- 
proaches the decay window, the worse is convergence of 
HH calculations. At energy E 2r = 360 kcV (which is al- 
ready close to three-body decay window E 3r = 344 keV) 
even FR procedure provides a convergence to the width 
value which is only about 65% of the exact value. 

In Fig. [5] the calculations with different E 2r values are 
summarized. The width grows rapidly as the two-body 
resonance moves closer to the decay window. The pen- 
etrability enhancement provided by the two-body reso- 
nance even before it moves into the three-body decay 
window is very important. Difference of widths with no 
core-p FSI and FSI providing the s-wave resonance to 
be at it experimental position E 2r = 0.535 MeV is more 
than two orders of the magnitude. The convergence of 
HH calculations also deteriorates as E 2r moves closer to 
the decay window. However, the disagreement between 
the HH width and the exact value is within the order 
of the magnitude, until the resonance achieves the range 
E 2r ~ (0.7 — 0.85)i?3 r . Within this range a transition 
from three-body to two-body regime happens (see also 
discussion in [8|), which can be seen as a drastic change 
of the width dependence on E 2r . This means that a se- 
quential decay via two-body resonance E 2r becomes more 
efficient than the three-body decay. In that case the hy- 
perspherical expansion can not treat the dynamics effi- 
ciently any more and the disagreement with exact result 
becomes as large as orders of the magnitude. The de- 
cay dynamics in the transition region is also discussed in 
details below. 

It can be seen in Fig. [5] that in three-body regime the 
dependence of the three-body width follows well the an- 
alytical expression 

T~ (E 3r /2-E 2r )- 2 (34) 

The reasons of such a behaviour will be clarified in 
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FIG. 9: Convergence of the 17 Ne width in a simplified model. 
Jacobi "T" system. Final state interaction describes s-wave 
p-p scattering. 



the forthcoming paper 37J ■ The deviations from this de- 
pendence can be found in the decay window (close to 
"transition regime") and at higher energies. This depen- 
dence is quite universal; e.g. for 45 Fe it is demonstrated 
in Fig. [TJJ where it follows the calculation results even 
with higher precision. 

Another important issue is a convergence of energy dis- 
tributions in the HH calculations, demonstrated in Fig. 
[6] for calculations with E 2r — 535 keV. The distribution 
is calculated in "Y" Jacobi subsystem, thus E x is the en- 
ergy between the core and one proton. The energy dis- 
tribution convergence is fast: the distribution is stable at 
K max = 10 — 14 and does not change visibly with further 
increase of the basis. There remain a visible disagree- 
ment with exact ("corrected") results, which give more 
narrow energy distribution. We think that this effect was 
understood in our work Q . The three-body calculations 
are typically done for p max ~ 500 — 2000 fm (p max ~ 1000 
fm everywhere in this work) . It was demonstrated in Rcf. 
[3] by construction of classical trajectories that we should 
expect a complete stabilization of the energy distribution 
in core+p subsystem at p max ~ 30000 — 50000 fm and the 
effect on the width of the energy distribution should be 
comparable to one observed in Fig. [6l 

The evolution of the energy distribution in core+p sub- 
system with variation of E 2r is shown in Fig. [7] When 
we decrease the energy E 2r the distribution is very stable 
until the two-body resonance enters the three-body de- 
cay energy window. After that the peak at about e ~ 0.5 
first drifts to higher energy and then for E 2r ~ 0.85E 3r 
the noticeable second narrow peak for sequential decay 
is formed. At E 2r ~ Q.7E 3r the sequential peak becomes 
so high that the three-body component of the spectrum 
is practically disappeared in the background. 

The result concerning the transition region obtained 
in this model is consistent with conclusion of the paper 
Q (where much simpler model was used for estimates). 
The three-body decay is a dominating decay mode, not 
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FIG. 10: Energy distributions for Ne in "T" Jacobi system 
(between two protons). 



only when the sequential decay is energy prohibited as 
Ei r > Ezr- Also the three-body approach is valid when 
the sequential decay is formally allowed (because E^r < 
E^r) but is not taking place in reality due to Coulomb 
suppression at E<i r > 0.8£ , 3 r . 

Geometric characters of potentials can play an impor- 
tant role in the width convergence. To test this aspect 
of the convergence we have also made the calculations 
for potential with repulsive core. This class of potentials 
was employed in studies of 17 Ne and 19 Mg in Ref. Q- A 
comparison of the convergence of HH calculations with s- 
wave 15 O+p potential from [tJ and Gaussian potential is 
given in Fig. [5] The width convergence in the case of the 
"complicated" potential with a repulsive core is drasti- 
cally worse than in the "easy" case of Gaussian potential. 
For typical dynamic calculations with K max = 20 — 24 the 
HH calculations provide only 20 — 25% of the width for 
potential with a repulsive core. On the other hand the 
calculations with both potentials provide practically the 
same widths T c [46| and FR provides practically the same 
and very well converged result in both cases. 



B. One final state interaction — p-p channel 

As far as two-proton decay is often interpreted as 
"diproton" decay we should also consider this case and 
study how important this channel could be. For this cal- 
culation we use a simple s-wave Gaussian p-p potential, 
providing a good low-energy p-p phase shifts, 



V(r) = -31exp[-(r/1.8) 2 



(35) 



Calculations with this potential are shown in Fig.[5](see 
also Table [Vj . First of all the penetrability enhancement 
provided by p-p FSI is much less than the enhancement 
provided by core-p FSI (the widths differs more than two 
orders of the magnitude, see Fig. [3]). This is the feature, 
which has been already outlined in our works. The p-p in- 
teraction may boost the penetrability strongly, but only 
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FIG. 11: Comparison of the OFSI calculations (solid lines) for 
45 Fe in the "T" system with diproton model Eq. (|36[) (dashed 
lines). Effective equivalent channel radius r c h(dp) for "dipro- 
ton emission" (a) as a function of radius po of the three-body 
potential (|25[) . the value po/y/2 should be comparable with 
typical nuclear sizes, (b) as a function of the position of the 
peak Ypoak in the three-body WF in Y coordinate. The 
dashe lines are given to guide the eye. 



in the situation, when protons occupy predominantly or- 
bitals with high orbital momenta. In such a situation the 
p-p interaction allows transitions to configurations with 
smaller orbital momenta in the subbarrier region, which 
provide a large increase of the penetrability. In our sim- 
ple model for 17 Ne 3/2~ state, we have already assumed 
the population of orbitals with minimal possible angular 
momenta and thus no strong effect of the p-p interaction 
is expected. 

Also a very slow convergence of the decay width should 
be noted in this case. For core-p interaction the K max ~ 
10 — 40 were sufficient to obtain a reasonable result. In 
the case of the p-p interaction the K max ~ 100 is required. 

Energy distributions between two protons obtained in 
this model are shown in Fig. 1101 Important feature of 
these distributions is a strong focusing of protons at small 
p-p energies. This feature is connected, however, not 
with attractive p-p FSI, but with dominating Coulomb 
repulsion in the core-p channel. This is demonstrated by 
the calculation with nuclear FSI turned off, which pro- 
vides practically the same energy distributions. Similarly 
to the case of the core-p FSI, very small K max > 10 is 
sufficient to provide the converged energy distribution. 
The converged HH distribution is very close to the exact 
("corrected") one but it is, again, somewhat broader. 

So far the diproton model has been treated by us as a 
reliable upper limit for three-body width Q . With some 
technical improvements this model was used for the two- 
proton widths calculations in Refs. 0, [l?], EH Ell- It 
is important therefore to try to understand qualitatively 
the reason of the small width values obtained in this form 
of OFSI model, which evidently represents appropriately 
formulated diproton model [47] . In Fig. [TT] we compared 
the results of the OFSI calculations for 45 Fe in the "T" 
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FIG. 12: Convergence of the 17 Ne width for experimental 
positions E 2r = 0.535 MeV of the + two-body resonance in 
the "X" subsystem and E 2r = 0.96 MeV of the 2+ two-body 
resonance in the "Y" subsystem (TFSI model). 



FIG. 13: Convergence of the Fe width for position of the 1 
two-body resonances in "X" and "Y" subsystems E 2r = 1-48 
MeV. 



Two final state interactions 



system with diproton width estimated by expression 



Td P = T7 2 i j \ Pi=o{0-$5E 3r ,r ch (dp) 7 2Z COIC ) , 

M TeA rl h {dp) 

(36) 

where M le d is the reduced mass for Ci-pp motion and 
fch(dp) is channel radius for diproton emission. The en- 
ergy for the relative 43 Cr-pp motion is taken 0.95E 3r bas- 
ing on the energy distribution in the p-p channel (see 
Fig. [TO] for example). In Fig. [TTk we show the effective 
equivalent channel radii for diproton emission obtained 
by fulfilling condition Td P = T c for OFSI model calcula- 
tions with different radii po 01 the three-body potential 
Eq. ([25]) . It is easy to see that for realistic values of these 
radii (pn ~ 6 fm for Fe) the equivalent diproton model 
radii should be very small (~ 1.5 fm). This happens pre- 
sumably because the "diproton" is too large to be con- 
sidered as emitted from nuclear surface of such small po 
radius. Technically it can be seen as the nonlinearity of 
the r c h(dp)-po dependence, with linear region achieved at 
pa ~ 15 — 20 fm. Only at such unrealistically large po val- 
ues the typical nuclear radius (when it becomes compa- 
rable with the "size" of the diproton) can be reasonably 
interpreted as the surface, off which the "diproton" is 
emitted. It is interesting to note that in the nonlinearity 
region for Fig. [TTk there exists practically exact corre- 
spondence between the Y coordinate of the WF peak in 
the internal region and the channel radius for diproton 
emission (Fig. [TTb). This fact is reasonable to interpret 
in such a way that the diproton is actually emitted not 
from nuclear surface (as it is presumed by the existing 
systematics of diproton calculations) but from the inte- 
rior region, where the WF is mostly concentrated. 



As we have already mentioned the situation of one final 
state interaction is comfortable for studies, but rarely re- 
alized in practice. An exception is the case of the El tran- 
sitions to continuum in the three-body systems, consid- 
ered in our previous work (23j . For narrow states in typ- 
ical nuclear system of the interest there are at least two 
comparable final state interactions (in the core-p chan- 
nel). For systems with heavy core this situation can be 
treated reasonably well as the Y coordinate (in "Y" Ja- 
cobi system) for such systems practically coincides with 
the core-p coordinate. Below we treat in this way 17 Ne 
(for which this approximation could be not very consis- 
tent) and 45 Fe (for which this approximation should be 
good). In the case of 17 Ne we are thus interested in the 
scale of the effect, rather in the precise width value. 

For calculations with two FSI for 17 Ne we used Gaus- 
sian d-wave potential (see Table [TJ , in addition to the 
s-wave potential used in Section IIV Al This potential 
provides a (i-wave state at 0.96 MeV (r = 13.5 keV), 
which corresponds to the experimental position of the 
first d-wave state in 16 F. The convergence of the 17 Ne 
decay width is shown in Fig. 1121 Comparing with Fig. 
[3] one can see that the absolute value of the width has 
changed significantly (2 — 3 times) but not extremely and 
the convergence is practically the same. Interesting new 
feature is a kind of the convergence curve "staggering" 
for odd and even values of K/2. Also the convergence 
of the corrected calculations requires now a considerable 

^max ^ 12 14. 

The improved experimental data for 2p decay of 45 Fe 
is published recently in Ref. pj]: E 3r = 1.154(16) MeV, 

r 2p = 2.85t°:6s x 10 ~ 19 MeV Pi/2(2p) = i.6±o:i H for 
two-proton branching ratio Br(2p) — 0.57. Below we use 
the resonance energy from this work. 

The convergence of the 45 Fe width is shown in Fig. [TO] 
The character of this convergence is very similar to that 
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FIG. 14: The 45 Fe g.s. width as a function of the two-body 
resonance position Ei r . Dashed, dotted and solid lines show 
cases of a pure HH calculation with K max = 24, the same but 
with Feshbach reduction from A" ma x = 100, and the corrected 
width r c . Gray area shows the transition region from three- 
body to two-body decay regime. The gray curve shows simple 
analytical dependence of Eq. (|34p . 
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FIG. 15: Interpolation of Ne decay width obtained in full 
three-body calculations by means of TFSI convergence curves 
(see Fig. [8}. Upper curves correspond to TFSI case with 
Gaussian potential in s-wave and compatible SI case for full 
three-body model. Lower curves correspond to TFSI case 
with repulsive core potential in s-wave and compatible GMZ 
case for full three-body model. 



in the 17 Ne case, except the "staggering" feature is more 
expressed. 

The dependence of the 45 Fe width on the two-body 
resonance energy E^ r is shown in Fig. [Mj Potential 
parameters for these 45 Fe calculations are given in Ta- 
ble HJ The result calculated for E 3r = 1.154 MeV and 
E% r = 1.48 MeV in paper 0] for pure [p 2 ] configuration 
is T = 2.85 x 1CT 19 MeV. The value K max = 20 was used 
in these calculations. If we take the HH width value from 
Fig. [Hat K max = 20 it provides T = 2.62 x 10" 19 MeV, 
which is in a good agreement with a full HH three-body 
model of Ref. (4J. However, from Fig. [13] we can conclude 
that in the calculations of 0] the width was about 35% 
underestimated. Thus the value of about T = 6.3 x 10~ 19 
MeV should be expected in these calculations. On the 
other hand much larger uncertainty could be inferred 
from Fig. Q3] due to uncertain energy of the 44 Mn ground 
state. If we assume a variation E^ r = 1.1 — 1.6 MeV the 



TABLE II: Parameters for 45 Fe calculations. Potential pa- 
rameters for p-wave interactions (|33[) in 43 Cr+p channel (V x o 
in MeV, r = 4.236 fm, r sph = 5.486 fm) and 44 Mn+p (V y0 in 
MeV, r = 4.268 fm, r sph = 5.527 fm), a = 0.65 fm. Calcula- 
tions are made with "effective Coulomb" of Eq. (|22[) . Widths 
T x , T y of the states in the subsystems are given in keV. Cor- 
rected three-body widths are given in the units 10 -19 MeV. 
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inferred from Fig. [TJ] uncertainty of the width would be 
T = (4 - 16) x 10" 19 MeV. On top of that we expect a 
strong p 2 1 f 2 configuration mixing which could easily re- 
duce the width within an order of the magnitude. Thus 
we can conclude that a better knowledge about spectrum 
of 44 Mn and a reliable structure information about 45 Fe 
are still required to make sufficiently precise calculations 
of the 45 Fe width. More detailed account of these issues 
is provided below. 

V. THREE-BODY CALCULATIONS 

Having in mind the experience of the convergence stud- 
ies we have performed large-basis calculations for Fe 
and 17 Ne. They are made with dynamical K max = 16—18 
(including Fechbach reduction from K max = 30 — 40) for 
17 Ne and K max = 22 (FR from K max = 40) for 45 Fe. 
The calculated width values are extrapolated using the 
convergence curves obtained in TFSI model (Figs. [15]) for 
17 Ne and [T3For 45 Fe). We have no proof that the width 
convergence in the realistic three-body case is asolutely 
the same as in the TFSI case. However, the TFSI model 
takes into account main dynamic features of the system 
causing a slow convergence, and we are expecting that 
the convergence should be nearly the same in both cases. 



A. Widths and correlations in Ne 

The potentials used in the realistic calculations are the 
same as used for 17 Ne studies in Refs. 0,111. The GPT 
potential (39j is used in the p-p channel. The core-p 
potentials are referred in [38[ as "GMZ" (potential in- 



13 




troduced in Q) and "high s" (with centroid of d-wave 
states is shifted upward which is providing a higher con- 
tent of s 2 components in the 17 Ne g.s. WF). Both poten- 
tials provide correct low-lying spectrum of 16 F and differ' 
only for d-wave continuum above 3 MeV (see Table IIIII) . 
The core-p nuclear potentials, including central, ss and 
Is terms, are taken as 



V(r) = 



V: 



(si • s 2 )V s < 



(!• 



2.0153 VL 



1 + cxp[(r - r l )/a\ 
x exp[(r - r l )/a] (l + exp[(r - r l Q )/a]y 2 ,(37) 

with parameters: a = 0.65 fm, rfj = 3.014 fm, rp >0 — 
2.94 fm, y c ° = -26.381 MeV, = -9 MeV, V c 2 = 
-57.6 (-51.48) MeV, V c 3 = -9 MeV, V s ° s = 0.885 MeV, 
V S 2 S = 4.5 (12.66) MeV, V ls = 4.4 (13.5) MeV (the values 
in brackets are for "high s" case). There are also repulsive 
cores for s- and p-waves described by a = 0.4 fm, rfj = 
0.89 fm, V covc = 200 MeV. These potentials are used 
together with Coulomb potential obtained for Gaussian 
charge distribution reproducing the charge radius of 15 O. 

To have extra confidence in the results, the width of 
the 17 Ne 3/2~ state is calculated in several models of 
growing complexity (Tables IIVHVIP . One can see from 
those Tables that improvements introduced on each step 
provide quite smooth transition from the very simple to 
the most sophisticated model. 

In Table HVl we demonstrate how the calculations in the 
simplified model of Section IIVI arc compared with calcu- 
lations of the full three-body model with appropriately 
truncated Hamiltonian. We can switch off correspond- 
ing interactions in the full model to make it consistent 
with approximations of the simplified model. To remind, 
the differences of the full model and simplified model are 
the following: (i) antisymmetrization between protons 
is missing in the simplified model and (ii) Y coordinate 
is only approximately equal to the coordinate between 
core and second proton. Despite these approximations 
the models demonstrate very close results: the worst dis- 
agreement is not more than 30%. 



In Table [V] we compare approximations of a different 
kind: those connected with choice of the Jacobi coordi- 
nate system in the simplified model. First we compare 
the "pure Coulomb" case: all pairwise nuclear interac- 
tions are off and the existence of the resonance is provided 
solely by the three-body potential (|25[) . This model pro- 
vides some hint what should be the width of the system 
without nuclear pairwise interactions. Then the models 
are compared with the nuclear FSIs added. The addition 
of nuclear FSI drastically increase width in all cases. It 
is the most "efficient" (in the sense of width increase) in 
the case of TFSI model in the "Y" system. Choice of 
this model provides the largest widths and can be used 
for the upper limit estimates. 

In Table IVT1 full three-body models are compared. The 
simplistic SI and S2 interactions correspond to calcula- 
tions with simplified spectra of the 16 F subsystem. For SI 
case it includes one s-wave state at 0.535 MeV (r = 18.8 
keV) and one d-wave state at 0.96 MeV (r = 3.5 keV). 
These are two lower s- and d-wave states known experi- 
mentally. In the S2 case we use instead the experimental 
positions of the higher component of the s- and d-wave 
doublets: s-wave at 0.72 MeV (r = 73.4 keV) and d- 



TABLE III: Low-lying states of 16 F obtained in the "GMZ" 
and "high s" core-p potentials. The potential is diagonal in 
the representation with definite total spin of core and proton 
S, which is given in the third column. 

Case GMZ high s Exp. 

r I S E 2r (MeV) T (keV) E 2r (MeV) T (keV) F (keV) 

r 

2" 
3" 
2" 
1" 









0.535 


18.8 


0.535 


18.8 


25(5) [34] 





1 


0.728 


73.4 


0.728 


73.4 


70(5) [34J 


2 





0.96 


3.5 


0.96 


3.5 


6(3) [34] 


2 


1 


1.2 


9.9 


1.2 


10.5 


< 15 [35] 


2 


1 


3.2 


430 


7.6 


~ 3000 




2 


1 


4.6 


1350 


~ 15 


~ 6000 
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nuclear pairwise potentials are turned off). 



wave at 1.2 MeV (r = 10 keV). Parameters of the core-p 
potentials can be found in Table [IJ Simple Gaussian p- 
p potential is used. The variation of the results 
between these models is moderate (~ 30%). The calcu- 
lations with GMZ potential provide the width for 17 Ne 
3/2 ~~ state which comfortably rests in between the re- 
sults obtained in the simplified SI and S2 models. The 
structure of the WF is also obtained quite close to these 
calculations. The structure in the "high s" case is ob- 
tained with a strong domination of the sd component. 
The width in the "high s" case is obtained somewhat 
larger (~ 11%) than in GMZ case, but this increase is 
consistent with the increase of the sd WF component, 
(~ 15%) which is expected to be more preferable for de- 
cay than d 2 component. 

It is important for us that the results obtained in the 
three-body models with considerably varying spectra of 
the two-body subsystems and different convergence sys- 
tematics appear to be quite close: r ~ (5 — 8) x 10~ 15 
MeV. Thus we have not found a factor which could lead 
to a considerable variation of the three-body width, given 
the ingredients of the model are reasonably realistic. 

The decomposition of the 17 Ne WF obtained with 
GMZ potential is provided in Table IVIII in terms of par- 
tial internal normalizations and partial widths. The cor- 
respondence between the components with large weights 
and large partial widths is typically good. However, there 
are several components giving large contribution to the 
width in spite of negligible presence in the interior. 

Complete correlation information for three-body de- 
cay of a resonant state can be described by two variables 
(with omission of spin degrees of freedom). We use the 
energy distribution parameter e = E x / E^ r and the angle 
cos((9fc) = (k. x ]ty) / (k x k y ) between the Jacobi momenta. 
The complete correlation information is provided in Fig. 
[T5]for realistic 17 Ne 3/2~ decay calculations. We can see 
that the profile of the energy distribution is characterized 
by formation of the double-hump structure, expected so 



far for p 2 configurations (see, e.g. Q)- This structure 
can be seen both in "T" system (in energy distribution) 
and in "Y" system (in angular distribution). In the cal- 
culations of ground states of the s-d shell nuclei we were 
getting such distributions to be quite smooth. It can be 
found that the profile of this distribution is defined by 
the sd/d 2 components ratio. For example in the calcula- 
tions with "high s" potential the total domination of the 
sd configuration leads to washing out of the double-hump 
profile. 

The correlations in the 17 Ne (shown in Fig. [TT5|) are 
strongly influenced by the nuclear FSIs. Calculations for 
only Coulomb pairwise FSIs left in the Hamiltonan are 



TABLE IV: Comparison of widths for 17 Ne (in 10~ 14 MeV 
units) obtained in simplified model in "Y" Jacobi system 
and in full three-body model with correspondingly truncated 
Hamiltonian. Structure information is provided for the three- 
body model. In the simplified model the weight of the [sd] 
configuration is 100% by construction. "No p-p" column 
shows the case where Coulomb interaction in p-p channel is 
switched off (see, (|2ip h "Eff." column corresponds to the ef- 
fective treatment (see, (|22|) ') of Coulomb interaction in the p-p 
channel in the simplified model, but to the exact treatment 
in full three-body model. 

pure Coulomb OFSI TFSI 





"no p-p" 


Eff. 


"no p-p" 


Eff. 


"no p-p" 


Eff. 


Simpl. 


0.017 


0.0032 


3.02 


0.545 


4.70 


1.37 


3-body 


0.024° 


0.0041° 


3.22 


0.555 


3.91 


0.445 


[sd] 


99.8 


99.3 


99.6 


99.5 


92.0 


72.6 


IP 2 } 


0.2 


0.6 


0.3 


0.4 


0.1 


0.2 


[d 2 ] 














7.8 


27.1 



"Small repulsion (~ 0.5 MeV) was added in that case in the p- 
wave core-p channel to split the states with sd and p 2 structure 
which appear practically degenerated and strongly mixed in this 
model. 
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(only s-wave core-p FSI). 



shown in Fig. 1171 The strong peak at small p-p energy 
is largely dissolved and the most prominent feature of 
the correlation density in that case is a rise of the distri- 
bution for cos(9k) — ► 1 in the "Y" Jacobi system. This 
kinematical region corresponds to motion of protons in 
the opposite directions from the core and is qualitatively 
understandable feature of the three-body Coulomb in- 
teraction (the p-p Coulomb interaction is minimal along 
such a trajectory). 

The distributions calculated in the simplified (OFSI) 
model are shown in Fig.[T5]on the same {e, cos(8k)} plane 
as in Figs. [TBI and [T71 It should be noted that here the 
calculations in "T" and "Y" Jacobi systems represent 
different calculations (with p-p FSI only and with core-p 
FSI only). In Figs. [TC] and [T7] two panels show differ- 
ent representations of the same result. Providing rea- 
sonable (within factor 2 — 4) approximation to the full 
three-body model in the sense of the decay width, the 
simplified model is very deficient in the sense of correla- 



TABLE V: Comparison of widths calculated for 17 Ne (1CT 14 
MeV units) and 45 Fe (10" 19 MeV units) with pure Coulomb 
FSIs and for nuclear plus Coulomb FSIs. Simplified OFSI 
model in "T", TFSI in "Y" Jacobi systems ("effective" 
Coulomb is used in both cases) and full three-body calcu- 
lations. 





pure Coulomb 


Nuclear+Coulomb 






"Y" 


3-body 




"Y" 


3-body 


17 Ne 


0.0011 


0.0032 


0.0041 


0.0077 


1.37 


0.76" 


M 


100 


100 


99.3 


100 


100 


73.1 


\P 2 } 








0.6 








1.8 


[d 2 ] 

















24.2 


45 Fe 


0.0053 


0.0167 


0.26 


0.034 


4.94 


6.3 6 



"This is a calculation with SI Hamiltonian. 
'This is a calculation providing pure p 2 structure. 



tions. The only feature of the realistic correlations which 
is even qualitatively correctly described in the simplified 
model is the energy distribution in the "Y" system. The 
"diproton" model (OFSI model with p-p interaction) fails 
especially strongly, which is certainly relevant to the very 
small width provided by this calculation. 



B. Width of 45 Fe 

The calculation strategy is the same as in [3] . We start 
with interactions in the core-p channel which give a res- 
onance in p-wave at fixed energy Ei r . Such a calculation 
provides 45 Fe with practically pure p 2 structure. Then 
we gradually increase the interaction in the /-wave, until 
it replaces the p-wave resonance at fixed i?2 r and then 
we gradually move the p-wave resonance to high energy. 
Thus we generate a set of WFs with different p 2 // 2 mix- 
ing ratios. 

The results of the improved calculations with the same 
settings as in |4( (the 44 Mn g.s. is fixed to have E^r = 1-48 
MeV) are shown in Fig. [19] (see also Table [Vj together 
with updated experimental data [l2j . The basis size used 
in [J] was sufficient to provide stable correlation pictures 
(as we have found in this work) and they are not updated. 



TABLE VI: Width (in 10" 14 MeV units) and structure of 
17 Ne 3/2~ state calculated in a full three-body model with 
different three-body Hamiltonians. 





SI 


S2 


GMZ 


high s 


J^max — 18 


0.35 


0.27 


0.14 


0.16 


Extrapolated 


0.76 


0.56 


0.69 


0.76 


[sd] 


73.1 


71.7 


80.2 


95.1 


[P 2 ] 


1.8 


1.8 


2.0 


1.3 


[d 2 ] 


24.2 


25.7 


16.8 


3.1 



1G 



The sensitivity of the obtained results to the experi- 
mentally unknown energy of Mn can be easily studied 
by means of Eq. (|34|) . The results are shown in Fig. [20] in 
terms of the regions consistent with experimental data 
on the {E 2r ,W(p 2 )} plane [IF(p 2 ) is the weight of p 2 
configuration in 45 Fe WF]. It is evident from this plot 
that our current experimental knowledge is not sufficient 
to draw definite conclusions. However, it is also clear 
that with increased precision of the lifetime and energy 
measurements for 45 Fe and the appearance of more de- 
tailed information on 44 Mn subsystem the restrictions on 
the theoretical models should become strong enough to 
provide the important structure information. 



VI. DISCUSSION 

General trends of the model calculations can be well 
understood from Tables IIVHVII For the pure Coulomb 
case the simplified model calculations (in the "Y" and 
"T" systems) and three-body calculations provide rea- 
sonably consistent results. The simplified calculations in 
the "Y" system always give larger widths than those in 
the "T" system. From decay dynamics point of view this 
leads to understanding of the contradictory fact that the 
sequential decay path is preferable even if no even virtual 
sequential decay is possible (as the nuclear interactions 
are totally absent in this case). 

The calculations with attractive nuclear FSIs rather 



TABLE VII: Partial widths IV7 of different components of 
17 Ne 3/2~ WF calculated in "T" Jacobi systems. Partial 
weights are given in "T" (value ) and in "Y" (value ) 
Jacobi systems. S x is the total spin of two protons. 



K 


L 




ly 




K'y 






2 


2 





2 





23.88 


33.87 


44.93 


2 


2 


2 








24.97 


16.52 


13.29 


2 


2 


1 


1 


1 


0.28 


7.39 


3.59 


2 


2 


1 


1 





1.54 






2 


2 





2 


1 


3.68 






2 


2 


2 





1 


3.68 






4 


2 





2 





8.97 


20.04 


3.19 


4 


2 


2 








8.68 


13.57 


5.57 


4 


2 


2 


2 





15.49 


0.32 


18.80 


4 


2 


1 


3 


1 


0.03 


2.18 


0.95 


4 


2 


3 


1 


1 





1.89 


0.63 


4 


1 


2 


2 


1 


1.02 






4 


2 





2 


1 


1.99 






4 


2 


2 





1 


2.07 






6 


2 


2 


4 





0.14 


0.77 


3.57 


6 


2 


4 


2 





0.14 


0.77 


0.78 


6 


2 





2 





0.50 


0.09 


0.69 


8 


2 


4 


4 





0.02 


0.003 


1.58 




£ 3r (MeV) 

FIG. 19: The lifetime of 45 Fe as a function of the 2p decay 
energy Es r . The plot is analogue of Fig. 6a from [4| with up- 
dated experimental data [12J and improved theoretical results. 
Solid curves shows the cases of practically pure p 2 and f 2 
configurations, dashed curves stand for different mixed p 2 // 2 
cases. The numerical labels on the curves show the weights 
of the s 2 and p 2 configurations in percents. 

expectedly provide larger widths than the corresponding 
calculations with Coulomb interaction only. The core- 
proton FSI is much more efficient for width enhancement 
than p-p FSI. This fact is correlated with the observation 
of the previous point and is a very simple and strong indi- 
cation that the wide-spread perception of the two-proton 
decay as "diproton" decay is to some extent misleading. 
As it has already been mentioned the p-p FSI influences 
the penetration strongly in the very special case when the 
decay occurs from high-Z orbitals (e.g. f 2 in the case of 
45 Fe). Thus we should consider as not fully consistent the 
attempts to explain two-proton decay results only by the 
FSI in the p-p channel (e.g. Ref. [19|]) as much stronger 
decay mechanism is neglected in these studies. 

From techical point of view the states considered in 
this work belong to the most complicated cases. The 
complication is due to the ratio between the decay energy 
and the strength of the Coulomb interaction (it defines 
the subbarrier penetration range to be considered dy- 
namically). Thus the convergence effects demonstrated 
in this work for 17 Ne have the strongest character among 
the systems studied in our previous works 0, H, 0, 
Because of the relatively small K max = 12 used in the 
previous works we have found an order of the magnitude 
underestimation of the 17 Ne(3/2 _ ) width. For systems 
like 48 Ni — 66 Kr the underestimation of widths in our 
previous calculations is expected to be about factor of 2. 
A much smaller effect is expected for lighter systems. 

It was demonstrated in [22], [23[ that the capture rate 
for the 15 0(2p,7) 17 Ne reaction depends strongly on the 
two-proton width of the first excited 3/2~ state in 17 Ne. 
This width was calculated in Ref. as 4.1 x 10~ 16 MeV 
(some confusion can be connected with misprint in Table 
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1.0 1.2 1.4 1.6 1.8 2.0 
E 2r (MeV) 

FIG. 20: Compatibility of the measured width of the 45 Fe 
with different assumptions about position i?2r of the ground 
state in the 44 Mn subsystem and structure of 45 Fe [weights 
of the p 2 configuration W(p 2 ) are shown on the vertical 
axis]. Central gray area corresponds to experimental width 
uncertainty T = 2.85±g;|| x 10" 19 MeV [H- The light 
gray area also takes into account the energy uncertainty 
E 3r = 1.154(16) MeV from [FJ]. The vertical dashed line 
corresponds to Ei r used in Fig. 1191 

III of Ref. 0|, see erratum). However, in the subsequent 
work [211, providing very similar to Q properties of the 
17 Ne WFs for the ground and the lowest excited states, 
the width of the 3/2" state was found to be 3.6 x 10" 12 
MeV. It was supposed in [2l[ that such a strong disagree- 
ment is connected with poor subbarrier convergence of 
the HH method in Q compared to Adiabatic Faddeev 
HH method of [2l|]. This point was further reiterated in 
Ref. [4l[. We can see now that this statement has a cer- 
tain ground. However, the convergence problems of the 
HH method are far insufficient to explain the huge dis- 
agreement: the width increase found in this work is only 
one order of magnitude. The most conservative upper 
limit T - 5 x 10" 14 MeV (see Table HVJ) was obtained in 
a TFSI calculation neglecting p-p Coulomb interaction. 
The other models systematically produce smaller values, 
with realistic calculations confined to the narrow range 
r - (5 - 8) x 1(T 15 MeV (Table Ell)- Thus the value 
r~4x 10" 12 MeV obtained in paper [21| is very likely 
to be erroneous. That result is possibly connected with a 
simplistic quasiclassical procedure for width calculations 
employed in this work. 

VII. CONCLUSION. 

In this work we derive the integral formula for the 
widths of the resonances decaying into the three-body 
channel for simplified Hamiltonians and discuss various 
aspects of its practical application. The basic idea of the 
derivation is not new, but for our specific purpose (pre- 



cision solution of the multichannel problem) several im- 
portant features of the scheme have not been discussed. 

We can draw the following conclusions from our stud- 
ies. 

(i) We presume that HH convergence in realistic calcu- 
lations should be largely the same as in the simplified 
calculations as they imitate the most important dynamic 
aspects of the realistic situation. The width values were 
somewhat underestimated in our previous calculations. 
The typical underestimation ranges from few percent to 
tens of percent for "simple" potential and from tens of 
percent to an order of magnitude in "complicated" cases 
(potentials with repulsive core). 

(ii) Convergence of the width calculations in the three- 
body HH model can be drastically improved by a simple 
adiabatic version of the Feshbach reduction procedure. 
For a sufficiently large dynamic sector of the basis the 
calculation with effective FR potential converges from 
below and practically up to the exact value of the width. 
For a small dynamic basis the FR calculation converges 
towards a width value smaller than the exact value, but 
still improves considerably the result. 

(iii) The energy distributions obtained in the HH calcu- 
lations are quite close to the exact ones. Convergence 
with respect to basis size is achieved at relatively small 
-Kmax values. The disagreement with exact distributions 
is not very significant and is likely to be connected not 
with basis size convergence but, with radial extent of the 
calculations 

(iv) Contributions of different decay mechanisms were 
evaluated in the simplified models. We have found that 
the "diproton" decay path is much less efficient than the 
"sequential" decay path. This is true even in the model 
calculations without nuclear FSIs (no specific dynamics), 
which means that the "sequential" decay path is some- 
how kinematically preferable. 

(v) The value of the width for 17 Ne 3/2" state was un- 
derestimated in our previous works by around an order of 
magnitude. A very conservative upper limit is obtained 
in this work as T ~ 5 x 10" 14 MeV, while typical values 
for realistic calculations are within the (5 — 8) x 10" 15 
MeV range. Thus the value r ~ 4 x 10" 12 MeV obtained 
in papers [U, |4l[ is likely to be erroneous. 

From this paper it is clear that the convergence issue 
is sufficiently serious, and in some cases were underesti- 
mated in our previous works. However, from practical 
point of view, the convergence issue is not a principle 
problem. For example the uncertain structure issues and 
subsystem properties impose typically much larger uncer- 
tainties for width values. For heavy two-proton emitters 
(e.g. 45 Fe) the positions of resonances in the subsystems 
are experimentally quite uncertain. For a moment this is 
the issue most limiting the precision of theoretical predic- 
tions. We have demonstrated that with increased preci- 
sion the experimental data impose strong restrictions on 
theoretical calculations allowing to extract an important 
structure information. 
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